Large-N transition temperature for superconducting films in a magnetic field 
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We consider tlie A'^-component Ginzburg-Landau model in the large A'^ limit, the system being 
embedded in an external constant magnetic field and confined between two parallel planes a 
distance L apart from one another. On physical grounds, this corresponds to a material in the 
form of a film in the presence of an external magnetic field. Using techniques from dimensional 
and zeta-function regularization, modified by the external field and the confinement conditions, we 
investigate the behavior of the system as a function of the film thickness L. This behavior suggests 
the existence of a minimal critical thickness below which superconductivity is suppressed. 
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PACS number(s): 74.20.-z, OS.lOCc, 11.25.Hf 



It is usually assumed that it is a good approximation to 
neglect magnetic thermal fluctuations in the Ginzburg- 
Landau (GL) model, when applied to study type II su- 
perconductors. This problem has been investigated by a 
number of authors, both in its single component and in its 
N- component versions. An account on the state of the 
subject can be found for instance in refs. P,|2iSSIEI3- 
In particular in ref.Q a large- iV theory of a second or- 
der transition for arbitrary dimension D is presented and 
the fixed point effective free energy describing the tran- 
sition is found. Here we investigate a confined version 
of the model studied in ref.Q- We consider the vector 
TV -component Ginzburg-Landau model in presence of an 
external magnetic field at leading order in , the system 
being submitted to the constraint of being confined be- 
tween two parallel planes a distance L apart from one 
another. Studies on confined field theory have been done 
in the literature since a long time ago. In particular, an 
analysis of the renormalization group in finite size ge- 
ometries can be found in rcf. 7] . This study is performed 
using a modified Matsubara formalism to take into ac- 
count boundary effects on scaling laws. From a physical 
point of view, for D = 3 and introducing temperature by 
means of the mass term in the Hamiltonian, the model 
studied here should correspond to a film-like material in 
presence of a magnetic field. We investigate the behavior 
of the system as a function of the separation L between 
the planes (the film thickness), using an extended com- 
pactification formalism in the framework of the effective 
potential, introduced in recent publications jl, ■ In 
this formalism has been employed to perform a study 
of the large- TV /3-function for superconducting films in 
a magnetic field. For the Ginzburg parameter k ^ 1 
(which is the case for high temperature superconductors) 
the Hamiltonian density of the GL model in an external 
magnetic field can be written in the form. 



H= |(V-zeA)0|2 



where V x A = H and 



(1) 



= a(T — To), with a > 
and To corresponding to the bulk transition temperature 



in absence of external field. This Hamiltonian density 
describes superconductors in the extreme type II limit. 
In the following we assume that the external magnetic 
field is parallel to the z axis and that the gauge A = 
(0, xH, 0) has been chosen. We will consider the model 
with N complex components and take the large iV 
limit at Nu fixed. If we consider the system in unlimited 
space, the field should be written in terms of the well 
known Landau level basis. 
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where Xi,p„,p^(r) are the Landau level cingcnfunctions 
given with energy eigenvalues i?/(|p|) = |pp + (2^ + l)w + 
TOq and ijj — eH is the so called cyclotron frequency. In 
the above equation p is a (Z) — 2)-dimensional vector. 

Now, let us consider the system confined between 
two parallel planes, normal to the z-axis, a distance L 
apart from one another and use Cartesian coordinates 
r = (z,z), where z is a {D — 3)-dimensional vector, with 
corresponding momenta k = (fcz,q), q being a (13 — 3) 
-dimensional vector in momenta space. In this case, the 
model is supposed to describe a superconducting mate- 
rial in the form of a film. Under these conditions the 
generating functional of the correlation functions is writ- 
ten as. 



Z = / V(j)*V(j)exp 



dz 



(3) 

with the field (l>{z,z) satisiying the condition of confine- 
ment along the z-axis, ip{z = 0,z) = ip(z ~ L.z). Then 
the field representation (0) should be modified and have 
a mixed series-integral Fourier expansion of the form. 
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where w„ = 2TTn/L, the label / refers to the Landau levels, 
and the coefficients c„ and 6(q) correspond respectively 
to the Fourier series representation over z and to the 
Fourier integral representation over the _D— 3-dimensional 
z-space. The above conditions of confinement of the z- 
dependence of the field to a segment of length L, allow us 
to proceed with respect to the z-coordinate, in a manner 
analogous as it is done in the imaginary-time Matsubara 
formalism in field theory. The Feynman rules should be 
modified following the prescription, 
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We emphasize that here we are considering an Euclidean 
field theory in D purely spatial dimensions, we are not 
working in the framework of finite temperature field the- 
ory. Temperature is introduced in the mass term of 
the Hamiltonian by means of the usual Ginzburg-Landau 
prescription. 

We consider in the following the zero external- 
momenta four-point function, which is the basic object 
for our definition of the renormalized coupling constant. 
The four-point function at leading order in is given by 
the sum of all chains of single one- loop diagrams. This 
sum gives for the L and lo -dependent four-point function 
at zero external momenta at the lowest Landau level ap- 
proximation the formal expression. 
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where !](£', L,lj) is the Feynman integral corresponding 
to the single one-loop subdiagram. 
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The sum over n and the integral over k can be treated 
using the formalism developed in 0| , which we resume 
below, adapted to the situation under study. The start- 
ing point is an expression of the form, 
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where we have used dimensionless quantities, — {m? + 
a;)/(47r^/i^) , [L^Y = ^ s-i^d is a mass scale. 
Note that our formalism makes sense only for dimensions 
D > 3. Using a well-known dimensional regularization 
formula I?] , Eq. (jSJ can be written in the form, 
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and Zf (s 
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; a) is one of the Epstein- Hurwitz zeta 



-functions defined by. 
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valid for Re{u) > K/2 (in our case i?e(s) > (D - 2)/2). 
The Epstein-Hurwitz zeta-function can be extended to 
the whole complex s-plane and we obtain, after some 
rather long but straightforward manipulations, 



U = h{D,s)uj 
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where 
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and are the Bessel functions of the third kind. 
Applying formula H12|l with s = 2 to the integral giving 
S(_D, L, cj), the result is that, we can write 2(1), uj) in 
the form, 

i:{D, L,uj)=uj [H{D, uj) + G{D, L, uj)] , (14) 

where the L and lo dependent contribution G{D, L,uj) 
comes from the second term between brackets in Eq. (|12|l , 
that is. 
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and H(D, w), is a polar parcel coming from the first term 
between brackets in Ea. lfT^ . 
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We see from Ea. (|16f) that for even dimension D ^ 6, 
H{D,uj) is divergent, due to the pole of the F-function. 
Accordingly this term must be subtracted to give the 
renormalized single bubble function S]jj(Z?, L, to). We get 
simply. 



^r{D,L,uj) =LuG{D,L,Lu). 



(17) 



For sake of uniformity the term H{D,uj) is also sub- 
tracted in the case of lower dimensions D, where no poles 
of F- functions are present. In these cases we perform a fi- 
nite renormalization. From the properties of Bessel func- 
tions, it can be seen from Ea. (|15|) that for any dimension 
D, G{D,L,uj) satisfies the conditions. 



lim G{D,L,u}) = 
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and G{D, L) > for any values of D, L. 

Taking inspiration from Eq. © , let us define the L and 
uj dependent renormalized coupling constant uu^D, L, lo) 
at the leading order in 1/N as, 

rg;^(0,L,c^) = -ur{D,L,oj) = 

—u 
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From Eas. (|19|l and (|17|l we can write the new L and w 
dependent renormalized coupling constant, 
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We see from Ea. (fTH|) that (3 = Nu corresponds to the 
renormalized coupling constant in absence of boundaries 
and of external field. Note that since the coupling con- 
stant (3 has mass dimension of 4 — D, then the coupling 
constant g{D,L,(jj) has mass dimension of 6 — £>. 

In order to study the critical behavior of our system, 
we start from the gap equation in absence of external 
field and of boundaries, in the disordered phase. 
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where ^ is the correlation length. In our case, in partic- 
ular in the neighborhood of the critical curve, the gap 
equation reduces to a (L,a;)-dependent Dyson-Schwinger 
equation. So, the generalization to our case of Eg. 1)21(1 
in the neighborhood of criticality can be written in the 
form, 
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In Ea. H22|l . we remember that = ■m'^{L,uj) + cu (the 
pole of the propagator of in presence of a magnetic 
field is at m^(i,cj) ~ — w)|[2, Ql and g{D,L,uj) is the 
renormalized (L, ti;)-dependent coupling constant, which 
is itself a function of via the mass m{L,uj). Per- 
forming steps analogous to those leading from Eq.© to 
Eq.lO, Eq.^ becomes 
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The couphng constant g{D,L,uj) is given by Ea. l(^ 
with G{D,L,Lu) in Ea. l(T5|l replaced by. 
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Eas. H23|l . (|20() and H24() are in fact a complicated set of 
integral equations involving since g{D,L,uj) is also 
dependent on the L and uj inverse correlation length. 
If we limit ourselves to the neighborhood of criticality, 
« 0, we may investigate the behavior of the system 
by using in both Eqs.lj^OI) and and in Eq.l(^ an 
asymptotic formula for small values of the argument of 
Bessel functions. 
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which allows after some straightforward manipulations, 
in the large TV -limit, to write Ea. H23l) in the form 
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where ({D — A) is the Riemann zeta-function, C(D — 4) = 
^^^(1/n^^'*), defined for D > 5 and we have used the 
label c to indicate that we are in the region of critical- 
ity). Similarly inserting Eq.(|5Sll in Eqs.(|2n|) and 
g{D,L,uj) can be written for « as. 
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We can not obtain critical lines in dimension D < 5 
by a limiting procedure from Ea. (|26|l . For D = 3, which 
corresponds to the physically interesting situation of the 
system confined between two parallel planes embedded in 
a 3-dimensional Euclidean space, we can obtain critical 
lines, by performing an analytic continuation of C(z) to 
values of the argument z < 1, by means of the reflection 
property of zeta-functions. 
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which defines a meromorphic function having only one 
simple pole at z = 1. We obtain taking mg = a{T — Tq), 
for £> = 3 (the physical dimension), the critical surface 

a(r, - To) +iu, + ^g{D = 3, L, luc)LC{2) = 0. (29) 

Using Eqs.(|2ZIl and |2H1) to evaluate g{D = 3,L,ujc), the 
tabulated values to the several F and C functions appear- 
ing in the above formulas, we obtain 
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a and (3 being the phenomenological Ginzburg-Landau 
parameters. Notice that as L cxd we obtain from Eq. 
the critical hne a{Tc - To) + ojc = for the bulk. 
In terms of the dimensionless quantities, respectively 
reduced critical field, temperature and film thickness, 
h = uJcil, * = Tc/Tq and I = L/^o, where ^0 = 
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(aTo)-i/2 is the zero-temperature Ginzburg-Landau co- 
herence length, the above equation can be rewritten as 

Hl^t) ^ \ -GOBI - Btl^ + Bl^ - 2880t: 

+ [{mBl + Btl^ - Bp + 28807r)2 
-I- 115207rB/3(l-t)]^/^|, (31) 

where B — /SS^q. The surface h — h{l,t) is illustrated 
in Fig. 1. We recall that, since we have used the low- 
est Landau level approximation in our calculations, this 
surface is only meaningful for high values of the external 
field, that is for low temperatures and thick films. 
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FIG. 1; Plot of the surface h = h{l, t) as defined by Eq. (glj, 
taking B = 10''. 

We can see from Fig. 1 that each value of I defines a 



critical line on the hy.t plane, corresponding to a film of 
thickness L. This set of critical lines suggests the exis- 
tence of a minimal value for the thickness L below which 
superconductivity is suppressed. Indeed, this can be seen 
from the plot of the reduced critical field at zero temper- 
ature, /iQ, as a function of the inverse of the reduced film 
thickness, shown in Fig. 2. This behavior may be con- 
trasted with the linear decreasing of Tc with the inverse 
of the film thickness in absence of external field that has 
been found experimentally in materials containing tran- 
sition metals, for example, in Nb and in W-Re alloys 
for these cases, it has been explained in terms of 
proximity, localization and Coulomb-interaction effects. 
Notice, however, that our results do not depend on mi- 
croscopic details of the material involved nor account for 
the influence of manufacturing aspects, like the kind of 
substrate on which the film is deposited. In other words, 
our results emerge solely as a topological effect of the 
compactification of the Ginzburg-Landau model in one 
direction, 
ho 
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FIG. 2: Zero-temperature reduced critical field as a function 
of the inverse of the reduced film thickness, obtained from Eq. 
lint taking f = and B = 10''. 
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